Operational quantification of continuous variable correlations 
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We quantify correlations (quantum and/or classical) between two continuous variable modes in terms of how 
many correlated bits can be extracted by measuring the sign of two local quadratures. On Gaussian states, such 
'bit quadrature correlations' majorize entanglement, reducing to an entanglement monotone for pure states. 
For non-Gaussian states, such as photonic Bell states, ideal and real de-Gaussified photon-subtracted states, and 
mixtures of pure Gaussian states, the bit correlations are shown to be a monotonic function of the negativity. This 
yields a feasible, operational way to quantitatively measure non-Gaussian entanglement in current experiments 
by means of direct homodyne detection, without a full tomographical reconstruction of the Wigner function. 
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Quantum information with continuous variables (CVs), re- 
lying on quadrature entanglement as a resource, has wit- 
nessed rapid and exciting progresses recently, also thanks to 
the high degree of experimental control achievable in the con- 
text of quantum optics |Q]]. While Gaussian states (coherent, 
squeezed, and thermal states) have been originally the pre- 
ferred resources for both theoretical and practical implementa- 
tions, a new frontier emerges with non-Gaussian states (Fock 
states, Schrodinger's cats, ...). The latter can be highly non- 
classical, possess in general more entanglement than Gaussian 
states [2], and are useful to overcome some limitations of the 
Gaussian framework such as entanglement distillation yfl and 
universal quantum computation [4]. Therefore, it is of central 
relevance to provide proper ways to quantify non-Gaussian 
entanglement in a way which is experimentally accessible. 

At a fundamental level, the difficulty in the investigation of 
entanglement - quantum correlations - can be traced back to 
the subtle task of distinguishing it from classical correlations 
iH. Correlations can be regarded as classical if they can be in- 
duced onto the subsystems solely by local operations and clas- 
sical communication, necessarily resulting in a mixed state. 
On the other hand, if a pure quantum state displays correla- 
tions between the subsystems, they are of genuinely quantum 
nature (entanglement). We adopt here a pragmatic approach: 
if two systems are in toto correlated, then this correlation has 
to be retrieved between the outcomes of some local measure- 
ments performed on them. We, therefore, investigate quadra- 
ture correlations in CV states. We are also motivated by the 
experimental adequacy: field quadratures can be efficiently 
measured by homodyne detection, without the need for com- 
plete state tomography. Specifically, in this paper, we study 
optimal correlations in bit strings obtained by digitalizing the 
outcomes of joint quadrature measurements on a two-mode 
CV system. First, we apply our procedure to Gaussian states 
(GS), finding that bit quadrature correlations provide a clear- 
cut quantification of the total correlations between the two 
modes. They are monotonic with the entanglement on pure 
states, and can be arbitrarily large on mixed states, the lat- 
ter possibly containing arbitrarily strong additional classical 



correlations. We then address non-Gaussian states (NGS), for 
which the exact detection of entanglement generally involves 
measurements of high-order moments [6]. The underlying 
idea is that for NGS obtained by de-Gaussifying an initial pure 
GS and/or by mixing it with a totally uncorrected state, our 
measure based entirely on second moments is still expected to 
be a (quantitative) witness of the quantum part of correlations 
only, i.e. entanglement. We show that this is indeed the case 
for relevant NGS including photon-subtracted states, photonic 
Bell states, and mixtures of Gaussian states. Notably, the com- 
plete entanglement picture in a recently demonstrated coher- 
ent single-photon-subtracted state J3l is precisely reproduced 
here in terms of quadrature correlations only. Our results ren- 
der non-Gaussian entanglement significantly more accessible 
in a direct, practical fashion. 

Quadrature measurements and bit correlations. — We con- 
sider a bipartite CV system of two bosonic modes, A 
and B, described by an infinite-dimensional Hilbert space. 
The probability distribution associated to the measurement 
of the rotated position quadrature Xa{0) in mode A with 
outcome x e A and uncertainty a is given by Va{%a) — 
Tr [p A R A (9)a(x A )R A (6y} = Tr [R A {0)^p A R A (e)a(x A )}, 
where &(xa) is a single-mode Gaussian (squeezed) state with 
first moments {x A , 0} and covariance matrix diagjcr 2 , 1/cr 2 }. 
Here Ra{&) is a unitary operator describing a rotation of 9 
on mode A, corresponding to a symplectic transformation 
given by R A (9) = ( c °f e ~^ n J) §]. Thus, one can either 
measure the rotated quadrature on the state (passive view) 
or antirotate the state and measure the unrotated quadrature 
(active view). Similarly we define Vb{^ V b) for mode B. 
The probability distribution associated to a joint measure- 
ment of the rotated quadratures Xa(0) and Xb(<p), is given 
by Vab(A^ b ) = Tr[p AB (R A (0) ® R B {tp))[t{x A ) ® 

a{x B )){R A {ey ®R B (^)\- 

We digitalize the obtained outputs by assigning the bits 
+ (— ) to the positive (negative) values of the measured 
quadrature. This digitalization transforms each joint quadra- 
ture measurement into a pair of classical bits. A string of 
such correlated bits can be used e.g. to distill a quantum 
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key J2I [Toll . Let us adopt a compact notation by denot- 
ing (at given angles 0, tp) = V A (±\x e A \), and P±| = 
(i I , =F \%b I) ■ The conditional probability that the bits 
of the corresponding two modes coincide is given by P AB = 
(Fab + ^ab)/E{ q =±,/3=±}^as- Correspondingly, the 
conditional probability that they differ is P AB = (V AB + 

Pi£)/E{a=±,/j=± } ^ Triviall y' p ab + pf B = i- if 

Pab > Pab tne measurement outcomes display correlations, 
otherwise they display anticorrelations. Notice that, if the two 
modes are completely uncorrected, P AB = P AB = 1/2. For 
convenience, we normalize the strength of bit correlations as 

B{\x e A \, \x%\) = 2\Px B - 1/2| = \PX B -Pf B \, (D 

so that for a completely uncorrected state E(|a;^J, \x B \) = 

0. The interpretation of Eq. (HJ in terms of correlations is 
meaningful if a fairness condition is satisfied: on each single 
mode, the marginal probabilities associated to the outcomes 
"+" or "-" must be the same: V\ = V A . V B = V E iflHl 

For a two-mode CV system, whose state is described by 
a Wigner function W, we define the 'bit quadrature corre- 
lations' Q as the average probability of obtaining a pair of 
classically correlated bits (in the limit of zero uncertainty) op- 
timized over all possible choices of local quadratures 

Q0S)=sup f f dx A dx B W(x A ,x B )[limB°(\x A \,\x B \)}, 

(2) 

where W{x A ,x B ) = ffdp A dp%W(x A7 p A ,x*,p%) is the 
marginal Wigner distribution of the (rotated) positions, and 
{x A ,p A ,x B ,p^} = {R(9) e R{(fi)){x A ,p A ,x B ,p B }. 
After some algebra, we can rewrite Eq. (f2]l as 
Q(p) = supe^E^ip)], where E%(p) = 
JJ dx A dx^sgn(x A x B )W(x A ,x B ) is the 'sign-binned' 
quadrature correlation function, which has been employed 
e.g. in proposed tests of Bell inequalities violation for CV 
systems 01211 . While this form is more suitable for an analytic 
evaluation on specific examples, the definition Eq. (|2]i is 
useful to prove the following general properties of Q: 

Lemma 1 (Normalization) : < Q(p) < 1. 
Proof. It follows from the definition of Q{p), as both B and 
the marginal Wigner distribution range between and 1. □ 
Lemma 2 (Zero on product states) : Q(p A (£> p B ) = 0. 
Proof. For a product state the probabilities factorize 

1. e. P AB = T" A P B and so P AB — P AB , where we have used 
the fairness condition. Namely B = 0, hence the integral in 
Eq. (O trivially vanishes. □ 
Lemma 3 (Local symplectic invariance) : Let U AtB be a uni- 
tary operator amounting to a single-mode symplectic opera- 
tion S AfB on the local phase space of mode A, B |8[]. Then 
QWa ® U B )p{(j\ ® U B )] = Q(p), 

Proof. Any single-mode symplectic operation S can be de- 
composed in terms of local rotations and local squeezings 
(Euler decomposition). By definition Eq. <J2J is invariant un- 
der local rotations, so we need to show that local squeez- 
ings, described by symplectic matrices of the form Z r = 



diag{l/r, r}, also leave Q(p) invariant. Adopting the pas- 
sive view, the action of local squeezings on the covariance 
matrix of each Gaussian state &(x A ' B ) is irrelevant, as we 
take eventually the limit a — > 0. The first moments are trans- 
formed as d A>B 1— ► Z~\d A ^ B , so that B AB ^(\x A \, \x B \) 1— > 
B A ~ B y °(\sx A \, \tx B \). On the other hand, the Wigner distribu- 
tion is transformed as W(g) h-> W^Z^ 1 Sum- 
ming up, local squeezings transform £ = {x Al p Al x B ,p B } 
into £ st = {sx A ,p A /s,tx B ,p B /t}. As Eq. ® involves in- 
tegration over the four phase space variables <i 4 £, we change 
variables noting that d 4 £ = d i £ st , to conclude the proof. □ 

It follows from Lemmas (1-2) that if Q > 0, then the state 
necessarily possesses correlations between the two modes. 
Lemma 3, moreover, suggests that Q embodies not only a 
qualitative criterion, but might be interpreted as bona fide op- 
erational quantifier of CV correlations. We will now show that 
this is the case for various important classes of states. 
Gaussian states. — Even though entanglement of GS is al- 
ready efficiently accessible via their covariance matrix |l|], we 
use such states as 'test-beds' for understanding the role of Q 
in discriminating CV correlations. The covariance matrix 7 of 
any two-mode GS p can be written in standard form as: 

= ( a 5 \ a = X a h, 13 = \ b h, 
T \5 T f3 )> 8 = diagK, -c p }, 

where, without loss of generality, we adopt the convention 
Cx > \cp\- The covariance matrix 7 describes a physical state 
if K.b > 1 and A < 1 + Det 7, with A = Det a + Det (3 + 
2Det7. The negativity 11411 . quantifying entanglement be- 
tween the two modes, reads N(p) — max{0, (1 — v)/(2v)}, 
where v" 1 = [A- (A 2 -4Det 7 ) 1 / 2 ]/2 with A = A-4Det5. 
For two-mode GS, Eq. (ffji evaluates to: 

Q(p) = (2 / 7r) arctan(c x / ^\ a \b - c 2 x ) , (3) 

where the optimal quadratures are the standard unrotated po- 
sitions (6 = tp = 0). First, we notice that Q = p 
describes a product state: for GS, Q > is then necessary 
and sufficient for the presence of correlations. Second, we 
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FIG. 1: (color online) Bit quadrature correlations vs normalized neg- 
ativity for 18000 random two-mode Gaussian states. The lowermost 
dotted (red) curve accommodates pure states. The leftmost solid 
(blue) vertical line denotes separable states parametrized by c p = 0, 
c x — e(A 2 — 1)/A, with < e < 1 and \ a ,t = A — > 00. The up- 
permost dashed (green) horizontal line denotes perfectly correlated 
states parametrized by c x = (A 2 — 1)/A, c p = ec x (and e, A as 
before). Product states (totally uncorrelated) lie at the origin. 
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observe that for pure GS [reducible, up to local unitary opera- 
tions, to the two-mode squeezed states p r = \<f> r )(<f> r \ charac- 
terized by A Qj b = cosh(2r) and c x = c p = sinh(2r)], Eq. 
yields a monotonic function of the negativity (see Fig. [TJ: Q 
is thus, as expected, an operational entanglement measure for 
pure two-mode GS. Third, we find that for mixed states Q ma- 
jorizes entanglement. Given a mixed GS />n with negativity 
N, it is straightforward to see that Q(pn)< [Eq. ®], is always 
greater than Q(\ip N )(ip N \) = (2 /it) arctan[(i> 2 - l)/(2z>)], 
with \ipjy) being a pure two-mode squeezed state with the 
same negativity N. Hence Q quantifies total correlations, and 
the difference Q{pn) — Q{4>n) (where the first term is due 
to total correlations and the second to quantum ones) can be 
naturally regarded as an operational measure of classical cor- 
relations II 1511 . We have evaluated Q on random two-mode 
GS as a function of their entanglement, conveniently scaled 
to 2N/(l + 2N), as shown in Fig. [TJ Note that for any en- 
tanglement content there exist maximally correlated GS with 
Q = 1, and also that separable mixed GS can achieve an ar- 
bitrary Q from to 1, their correlations being only classical. 

Non-Gaussian states. — Let us now turn our attention towards 
NGS, whose entanglement and correlations are, in general, 
encoded in higher moments as well. We focus on the most 
relevant NGS recently discussed in the literature and some- 
times experimentally realized. Remarkably, we find for all of 
them a monotonic functional dependence between the entan- 
glement (negativity) and the quadrature correlations Q. 
Photonic Bell states. We consider Bell-like states of the 
form |$±) = VPlOO) ± y/T=p\ll) and |*±) = ^01) ± 
\J\ — p|10) (with < p < 1), which are nontrivial exam- 
ples of superpositions of Fock states, entangled with respect 
to the (discrete) photon number. The negativity of these states 
is Nb — — p). For the four of them, Eq. (fJJ reads 

Qb = (4:/tt)Nb, showing a perfect agreement between the 
quadrature CV correlations and the entanglement. 
Photon subtracted states. Most attention is being drawn by 
those CV states obtained from Gaussian states via subtract- 
ing photons 10, EH]: they perform better as resources for 
protocols like teleportation IU7[ Hal and allow for loophole- 
free tests of nonlocality lfl3ll . Let us recall their preparation, 
following 11711 . The beam A (B) of a two-mode squeezed 



state \cf) r ) is let to interfere, via a beam splitter of transmit- 
tivity T, with a vacuum mode A' (£?')■ The output is a four- 
mode Gaussian state of modes AA'BB'. Detection of one 
photon in each of the two beams A' and B', conditionally 
projects the state of modes AB into a pure NGS, given in the 
Fock basis by \i> P V) = J2n=o C «K n )< where c„(A,T) = 
(n + 1)(TA)" (1 - T 2 A 2 ) 3/2 /VT 2 A 2 + 1, and A = tanh r. 
The limiting case T — *■ 1, occurring with asymptotically van- 
ishing probability, corresponds to an ideal two-photon sub- 
straction, \ip p V) oc (ciaclb)^)- For any T, the negativity is 
computable as N(ip p V) = 2/(1 - TA) - 1/(1 + T 2 A 2 ) — 1. 
It increases with r and with T: only in the case T = 1 it ex- 
ceeds N(<f) r ) for any r, diverging for r — > oo. For all T < 1, 
the entanglement of \ipp S ) eventually saturates, and above a 
squeezing threshold rr ^> 1, the original GS is more entan- 
gled than the resulting non-Gaussian one. Following B13I1 . the 
explicit expression of the quadrature bit correlations Eq. (0 
can be analytically obtained: 



Q(4V) = 



oo 

E 

n>m— 



2 m+n+3 TT[J r (m, n) - T{n, m)] 2 c m c n 



(m 



where the c fc 's are defined above, T = [V (-f ) T (^)] , 
and r is the Euler function. As depicted in Fig. |2j the be- 
haviour of the entanglement is fully reproduced by Q(ip ps ) 
which again is a monotonic function of N: Eq. (01 is thus 
measuring truly quantum correlations of this important class 
of NGS as well. 

Experimental de-Gaussified states. Up to now we considered 
nearly-ideal non-Gaussian situations. We apply now our def- 
inition to the coherent photon-subtracted state p e ™ recently 
studied and demonstrated by Ourjoumtsev et al. [7]. Refer- 
ring to their paper for details on the state preparation, we just 
mention that each of the two beams of the two-mode squeezed 
state \<f) r ) is let to interfere with the vacuum at a beam splitter 
with reflectivity R <C 1, and by using an avalanche photodi- 
ode a single photon is subtracted from the state in a delocal- 
ized fashion. The realistic description of the obtained highly 
mixed state involves many parameters: we fix all to the values 
obtained from the theoretical calculations and/or experimen- 
tal data of Ref. [7], but for the reflectivity R and the initial 
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FIG. 2: (color online) (a) Bit quadrature correlations Q vs squeez- 
ing r and beam splitter transmittivity T for photon-subtracted states 
(shaded surface) and two-mode squeezed states (wireframe surface), 
(b) Q vs normalized negativity for photon- subtracted states. 




FIG. 3: (color online) Bit quadrature correlations Q vs squeezing 
r (in decibels) for the de-Gaussified states demonstrated in [7J] with 
beam splitter reflectivities R equal to (from top to bottom): 3% (red), 
5% (blue), 10% (magenta), 20% (green). The dashed black curve 
depicts Q for a two-mode squeezed state. Cfr. Fig. 6 of [7J]. 
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squeezing r which are kept free. We then evaluate Eq. (O 
as a function of r for different values of R. Unlike the previ- 
ous cases, optimal correlations in the state p exp occur between 
momentum operators (8 = tp = tt/2). Also for this realistic 
mixed case, the correlation measure Q reproduces precisely 
the behaviour of the negativity, as obtained in |01 after full 
Wigner tomography of the produced state p exp . In particular, 
the negativity (and Q) increases with the squeezing r, and de- 
creases with R. Below a threshold squeezing which ranges 
around ~ 3 dB, the NGS exhibits more entanglement (larger 
Q) than the original two-mode squeezed state. Our results de- 
picted in Fig. [3] compare extremely well to the experimental 
results (Fig. 6 of [7]) where the negativity is plotted as a func- 
tion of r for different i?'s. 

Mixtures of Gaussian states etc. Recent papers lfl9lEoll dealt 
with mixed NGS of the form p m = p\ 4> r ){4> r | + ( 1 -p) | 00)(00 1 , 
with < p < 1. They have a positive Wigner function 
yet they are NGS (but for the trivial instances p — 0, 1). 
Clearly, the de-Gaussification here reduces entanglement and 
correlations in general. The negativity of such states reads 
N m {pm) = pN(4> r ) = p(e 2r — l)/2 and is increasing both 
with r and with p. The same dependency holds for the bit 

correlations, Q(p m ) = (2p/7r) arctan |iV„ 
which again is a monotonic function of N r7 



2N m +p 1 

for any p. We fur- 
ther studied other NGS including photon-added and squeezed 
Bell-like states lfl8ll . and their mixtures with the vacuum: for 
all we found a direct match between entanglement and Q i2lll . 
Discussion. — By analyzing the maximal number of corre- 
lated bits (Q) that can be extracted from a CV state via quadra- 
ture measurements, we have provided an operational quantifi- 
cation of the entanglement content of several relevant NGS 
(including the useful photon-subtracted states). Crucially, one 
can experimentally measure Q by direct homodyne detections 
(of the quadratures displaying optimal correlations only), in 
contrast to the much more demanding full tomographical state 
reconstruction. One can then easily invert the (analytic or nu- 
meric) monotonic relation between Q and the negativity to 
achieve a direct entanglement quantification from the mea- 
sured data. Our analysis demonstrates the rather surprising 
feature that entanglement in the considered NGS can thus be 
detected and experimentally quantified with the same com- 
plexity as if dealing with GS. 

Interestingly, this is not true for all CV states. By defini- 
tion, Q quantifies correlations encoded in the second canoni- 
cal moments only. We have realized that there exist also states 
[e.g. the photonic qutrit state \ip h ) = \00)/y/2 + (|02) + 
|20))/2] which, though being totally uncorrelated up to the 
second moments (Q = 0), are strongly entangled, with cor- 
relations embedded only in higher moments. The character- 
ization of such states is an intriguing topic for further study 
1 2 111 . In this respect, it is even more striking that the measure 
considered in this paper, based on (and accessible in terms 
of) second moments and homodyne detections only, provides 
such an exact quantification of entanglement in a broad class 
of pure and mixed NGS, whose quantum correlations are en- 



coded nontrivially in higher moments too, and currently repre- 
sent the preferred resources in CV quantum information. We 
focused on optical realizations of CV systems, but our frame- 
work equally applies to collective spin components of atomic 
ensembles [22], and radial modes of trapped ions 1I23I1 . 
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